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I. INTRODUCTION 



section. 



Stimulated by the developments of superstring the- 
ory, various cosmological solutions based on string theory 
have been proposed. Although the theory has not devel- 
oped enough to present a unique history of our universe 
in its earliest stage, one may now imagine that the big 
bang is no longer just a point of enigma named the initial 
singularity but has a rich and complex structure. 

Among these string-based universe models, the most 
widely studied would be the so-called pre-big-bang model 
HI (in the literature it is often called string cosmology). 
The remarkable aspect of this model is that it tries to 
explain the inflationary behavior of the early universe by 
introducing a polelike acceleration phase (superinflation) 
which is driven by the kinetic term of the dilaton. This 
superinflationary branch of the solution has a dual rela- 
tion named the scale factor duality with the Friedmann 
branch, which is the usual decelerating expansion of the 
universe. The biggest problem arising in the pre-big-bang 
model is the difficulty of connecting the superinflationary 
branch and the Friedmann branch (termed the graceful 
exit problem), and there are no-go theorems proved un- 
der some assumptions [||. 

Antoniadis, Rizos, and Tamvakis Q proposed a non- 
singular (i.e. free of the graceful exit problem) cosmo- 
logical model by including the 1-loop (genus) correction 
term in the low-energy string effective action. This 1- 
loop model is excellent in that it gives a simple exam- 
ple of the smooth transition between the superinflation- 
ary branch and the Friedmann branch. This original 
1-loop model involves dilaton and modulus fields, and 
there is a simplified version [|| by Rizos and Tamvakis 
with a modulus field only. This metric-modulus sys- 
tem gives essentially the same nonsingular solution as 
the full metric-dilaton-modulus system, since the behav- 
ior of the solution mainly depends only on the modulus 
field. The existence of the nonsingular solution is analyti- 
cally shown by Rizos and Tamvakis j| for the Friedmann- 
Robertson- Walker metric, i.e., assuming the homogeneity 
and isotropy of the universe. We studied whether the na- 
ture of this nonsingular solution is affected if anisotropy is 
included. The purpose of this paper is to extend the solu- 
tion of this metric-modulus system to include anisotropy 
and to observe the behavior of its solutions, particularly 
of nonsingular ones. 

The following sections of this paper are organized as 
follows. In the Sec. II we briefly review the isotropic 
case studied by Rizos and Tamvakis Q, and then derive 
the basic equations of motion from the effective action 
for the Bianchi type I metric. We also solve these equa- 
tions analytically in the asymptotic region. In Sec. Ill 
we solve these equations numerically. We study the so- 
lutions through several cross sections of the parameter 
space. The existence of nonsingular anisotropic solutions 
is shown, and the nature of the singularity is also exam- 
ined. Implications of our results are discussed in the last 



II. MODEL AND EQUATIONS OF MOTION 

We start with the action]] given by [Q 

S = j d^x^—g^R-^Daf- A^^j , (1) 

which is essentially the same as the 1-loop-corrected 4- 
dimensional effective action of orbifold-compactified het- 
erotic string ||, 

S = J d*xV=-g{lR-\mf -\{D*Y 

+ A[A ie *-A 2 C(a)]i?^|, (2) 

except that the dilaton field $ is neglected. R, <E>, and 
a are the Ricci scalar curvature, the dilaton, and the 
modulus field, respectively. Our convention is g^ v = 

{—, +, +, +), R^avf} = ^^afj.v + " "j Raj3 = R^a^P, and 

87rG = 1. The Gauss-Bonnet curvature is defined as 
r gb = R^^R^KX-iR^R^u + R 2 , and f(o-) is a func- 
tion determining the coupling of a and the geometry, 
written in terms of the Dedekind r\ function as 

£(a) = -ln[2eV(^)] 

tj oo 

= -ln2-a+^--4^1n(l-e- 2 "^). (3) 

3 n=l 

This ^(cr) , an even function of a, has a global minimum at 
cr = and increases exponentially as a — > ±oo. Ai is the 
4-dimensional string coupling and takes a positive value. 
A2 is proportional to the 4-dimensional trace anomaly of 
the N = 2 sector and determined by the number of chiral, 
vector, and spin-| supermultiplets. It is important that 
A2 can take positive values, since nonsingular solutions 
arise only when A2 > 0. In our simplified model (Q), 
therefore, we assume A to be positive (in actual numerical 
calculations we set A = 1) and adopt the form of the £ 
function (||). 

A. Isotropic solutions 

First, we review the homogeneous and isotropic case 
which is discussed in j| . We neglect the spatial curvature 
and write the metric in the flat FRW form 



1 This action is also justifiable in more solid grounds. In the 
N = 4 superstring compactifications (heterotic on T 6 or type 
II A on K3 x T 2 ) the function £(<r) given in Eq.(|^) describes 
the exact R 2 couplings including nonperturbative effects Q. 
These vacua exhibit an exact 5* duality, which is identical with 
the SL(2, Z) modular invariance on a. 
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ds 2 



-N{t) 2 dt 2 + a{t) 2 (dx 2 + dy 2 + dz 2 ). 



(4) 



Variation of the action ([[]) with respect to the lapse N, 
the scale factor a, and the modulus field a gives three 
equations of motion as 



i2// ->// 2 ; (1 ^//i 



= 0, 



a + 3H& 



3A 
T 



oa 



0, 



(5) 
(6) 
(7) 



where 7? is the Hubble parameter a/a, the overdot means 
derivative with respect to physical time t, and we have 
set N(t) = 1. As a result of the absence of scales [we are 
only considering spatially flat metric (Q)] and the exis- 
tence of the constraint (|j), the solutions are completely 
determined by a couple of first order differential equa- 
tions for two variables H and cr; that is, if we give values 
of H and a at some time t, the preceding and following 
evolutions of the solution are automatically determined 
by these equations. 

Figure 1 shows the H-a phase diagram of the isotropic 
system solved with initial conditions H > and a > 0. 
These solution flows are distinguished by only one degree 
of freedom (for example, the value of H at some fixed 
cr > 0). There are singular solutions and nonsingular 
solutions, and it is shown in Q that all flows in the H > 
0, a < quarter-plane continue smoothly to H > 0, 
a > quarter-plane, but some flows in H > 0, a > 
quarter-plane go into singularity and do not continue to 
the a < region. It is also shown in Q that the signs 
of H and a are conserved throughout the evolution of 
the system. Since a is always positive in Fig. 1, time 
flows from left to right. We consider that our Friedmann 
universe corresponds to the "future" region (cr > 0) in 
Fig. 1, and we regard the "past" region (cr < 0) with 
increasing Hubble parameter as a superinflation, which 
is expected to solve the shortcomings of the big-bang 
model. 



B. Equations of motion 

We now extend the above model to anisotropic Bianchi 
type I space-time. We write the metric as 

ds 2 = -N{t) 2 dt 2 + e 2a ^dx 2 + e 2 ^dy 2 + e 2 ^dz 2 , 

(8) 

and define the anisotropic expansion rates as 

p = a, q=P, r = 7. (9) 

The average expansion rate, which coincides with the 
Hubble parameter H in the isotropic limit, is 



#avr = ^(p + Q + r). 



(10) 



The equations of motion are obtained by variation of the 
action (hi) with respect to N, a, (3, 7, and cr, viz., 



(11) 
(12) 
(13) 
(14) 
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1 <9£ 

-A— {pqr + pqr + pqr + pqr(p + q + r)} = 0, (15) 



where 
A = l 



A ^,„ , A ' f^Y 2 



A <9£ . A 2 (dtV a 

n A ' f^Y 2 2 
A 2 {d{\\ 2 2 



E = — — — I r z p 

P A " W 2 2 



2 2<9cr 
2 



cr 



A ^„„„ A /^\ 22 



1.2 , A9 2 ^. 



rr 2 2 x • 2 , n w S . 2 

// = -«? -r -?r--<7 +2^r 
A ^_ x fdt\ 2 „2 2, 



and 



(16) 
(17) 
(18) 
(19) 
(20) 
(21) 



--— apqr- - I — ) p 2 q 2 r(p + q + r), (22) 



2 ~do° Pqr ~ 4 \da) M T ( p + q + r ^ ( 23 ) 

„ 2 2 1 • 2 A . n 

Q = -r -p -rp--a + - -^a rp 

A ^,„„„ x fdtY 2 2, 



2 ~b~a° Pqr ~ 4 I 9c7 ' P 1 r (p + 1 + r )> ( 24 ) 



A = 2 ABC + DEF - FC 2 - DA 2 - EB 2 . (25) 

In the isotropic limit, Eqs. ( pT| ) an d (p"5| ) become Eqs. 
d) and (0) respectively, and Eqs. ((Tg^ dTJ), and @) 
reduce to one equation (||). Similar to the isotropic case, 
once the values of cr, p, q, and r are given at some time 



3 



t, and a is given according to the constraint (|1 1\) , the 
system evolves along the flow of the solution following 
Eqs. (Ill) — (|l5|) . In our numerical analysis we integrated 
Eqs. (|12|) — (]15| ) with the initial condition satisfying the 
constraint (|l 1[) . This constraint equation (11) is also used 
to estimate the numerical error. 

There are five parameters (a, a, p, q, and r) and one 
constraint (0); so the solution flows are drawn in a 4- 
dimensional space, for example, a-p-q-r space. Since the 
aim of this paper is to examine the effect of anisotropy 
on the nonsingular solution, we observe the change of the 
solution as we increase the anisotropy of the metric. To 
facilitate this we introduce two parameters indicating the 
anisotropy of the metric, viz., 



A : p — — + uj 2 + ^3 = sgn(i), 



X = 



P 



Y = 



P- 



q — r 
p + q + r 



(26) 



In this notation, (X, Y) = (0, 0) corresponds to the 
isotropic (FRW) metric, and postulating i? avr to be pos- 
itive, the region surrounded by the lines Y = 2X + 1, 

Y = \X —\,Y = —X + 1 indicates the universe expand- 
ing in all directions (these lines are drawn in Figs. 3 and 
6). The regions Y < 2X + 1, Y > \X - §, Y > -X + 1, 
etc., are the "Kasner-like" universe expanding in two di- 
rections and contracting in one direction. The regions 

Y > 2X + 1, Y > \X-\, Y > -X + l, etc., describe the 
universe expanding in one direction, contracting in two 
directions. The universe shrinking in all directions can- 
not be included if we use Eqs. ( p6] ) and assume H avl > 0. 
Instead of a-p-q-r we use a-H avr -X-Y as four variables 
describing our system. 



C. Asymptotic solutions 

Before going to the numerical analysis, we study the 
asymptotic form of the solutions at t — > ±oo, \a\ — ► oo. 
We assume \a\ to become large when t — > ±oo. The 
asymptotic form of the derivatives of the function £ is 



da 
d 2 lj 
da 2 



sgn(cr)-e 



r 1 ' 



(27) 
(28) 



If we assume the power-law ansatz for the expansion 
rates, the asymptotic form of the modulus has to be loga- 
rithmic in order to cancel the exponential dependence of 
Eqs. (|27]) and (pq). Thus we choose the following forms 
for the asymptotic solutions: 



p~wi|t|", 
q~LJ 2 \t\ p , 
r~uj 3 \t\?, 
a ~ (To + k/4ln \t\ 



Putting all these into Eqs. (Ill 
possible asymptotic solutions 



(29) 
(30) 

(31) 
(32) 

151), we obtain two 



2 i 2 
UJ3 + <^4 



B:p = -2,|w 4 | = 5, 
^i^2^3 = -sgn(t) 



5exp[-o- sgn(cj 4 )] 
Att 



(33) 
(34) 



The solution A is obtained by balancing terms that do 
not include the 1-loop effect, i.e., the Gauss-Bonnet term, 
and thus describes the asymptotic behavior where the 
Gauss-Bonnet effect is negligible. In the absence of the 
modulus field, A is nothing but the Bianchi type I vac- 
uum (Kasner) solution. The solution B is obtained by 
balancing the kinetic term of the modulus and the Gauss- 
Bonnet term, and this solution corresponds to the phase 
where the Gauss-Bonnet term is important. Other pos- 
sibilities of solutions are excluded as long as we impose 
A > 0, which is, in the isotropic limit, a necessary condi- 
tion for the existence of the nonsingular solutions. 

Following the isotropic case , we choose a > and 
assume the solution A in the future asymptotic region 
and B in the past asymptotic region. Then, in the region 
t — » oo, the conditions ([||) can be seen in the 
space as the cross section of the sphere of radius y/l — w| 
centered at the origin with the u>i + io 2 + ^3 = 1 plane. 
Depending on the asymptotic value of W4, the asymptotic 
solutions of t — > 00 are categorized into two cases. 

Al: < uii < yf^. The asymptotic solution is either 

expanding in all directions (Friedmann type) or expand- 
ing in two directions, shrinking in one direction (Kasner 
type). 

A2: \J\ < 0J4 < y^|. The asymptotic solution is 
Friedmann type only. 

In terms of X and Y introduced in Eq. (|26|), the 
asymptotic solution is represented by a point on an arc of 
the ellipse X 2 + Y 2 - XY = 1 - |cj|. Thus A2 falls into 
the region inside the oval X 2 + Y 2 - XY = |, and Al 
is in the region between two ellipses X 2 + Y 2 — XY = \ 
and X 2 + Y 2 — XY = 1. No asymptotic solutions exist in 
the region outside the ellipse X 2 + Y 2 — XY = 1, where 
the constraint equation (|ll| ) is not satisfied. Therefore, 
at least in the far enough future region, it is sufficient 
to examine solutions near the isotropic one. The dif- 
ference in our model from the Kasner (Bianchi type I 
vacuum) solution is the existence of the field a, which 
allows the existence of a Friedmann-type solution in the 
future asymptotic region. 

In the past asymptotic region, the condition (|34|) indi- 
cates W1W2W3 > 0, i.e., one of the following. 

Bl: u>i, u> 2 , u>3 > 0. 

B2: One of u>i (i — 1, 2, 3) is positive; two are negative. 

This means that in the past asymptotic region the uni- 
verse is either expanding in all directions or expanding 
in one direction, contracting in two directions. 
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III. NUMERICAL RESULTS 

Once the anisotropy is included, the behavior of the 
solution deviates substantially from the isotropic case. 
Figure 2a shows the average expansion rate H avr = (p + 
<7 + r)/3 versus the modulus a in the anisotropic (Bianchi 
type I) case, solved with initial anisotropy X = 0.1, 
Y = 0.2 at a = —10. Unlike the isotropic case (Fig. 1), 
some solution flows in the a < region do not continue 
smoothly to the a > region, but terminate suddenly 
with finite values of a and H avv . At these unusual singu- 
larities the time derivatives of p, q, and r become infinite, 
although p, q, and r themselves stay finite (see Fig. 2b). 
This is because the value of A, Eg. ( p5|), approaches zero, 
while the numerators of Eqs. ([l2|)' (|13]), and ( |l4| ) stay 
finite. Thus, the function A plays an important role in 
the anisotropic case, and the regularity of the solutions 
depends largely on its behavior. 

In the equations in our model there are four indepen- 
dent variables a-H avl -X-Y . Since our interest is mainly 
in the vicinity of the isotropic solution, we examine the 
solutions which pass near the origin of the X-Y plane, 
first at the a = —10 section with several different values 
of H avl and next at the a = 2.5 section. 



A. Solutions through the a = —10 cross section 

It is helpful to consider the general behavior of A 
and the region prohibited by the constraint equation be- 
fore solving the equations numerically. Rewriting A, Eq. 
(p5|), and the constraint equation jTl| ) in terms of cr, H avr , 
X, and Y, we can specify A > 0, A < 0, and prohibited 
regions on the X-Y plane for fixed a and H avl , which is 
shown in Fig. 3. The dark-shaded region is the prohib- 
ited region, the light-shaded region is where A < 0, and 
the white region is where A > 0. Since A = is not 
physically allowed [p, q, and f diverge from Eqs 



(12) 



(|13|), and (|1J)], the solutions in the white region cannot 
go smoothly to the light-shaded region. Also indicated in 
Fig. 3 are the lines X+Y = 1, etc., discussed in Sec. II B. 
We can see that the universe expanding in all directions 
always lies in the A > region. For larger H avY the pro- 
hibited region becomes thinner, and the white and light 
gray regions will be separated by the lines X + Y = 1, 
etc. 

Starting from the initial conditions u = —10 and 
#avr = 0.001, 0.005, 0.01, we solved the equations fu- 
tureward and indicated the behavior of the solutions in 
the X-Y plane (Fig. 4) . Because of the symmetry of the 
axes, we restrict the region to X > 0, Y > 0, and since 
we are not interested in the prohibited region, we only 
examined the vicinity of the origin. 

The black region in Fig. 4 is prohibited by the Hamilto- 
nian constraint, and the regions marked NS means non- 
singular solutions. The difference between NSa and NSb 
is in their form in the future asymptotic region, where 



NSa has Friedmann-type (expanding in all directions) 
and NSb has Kasner-type (expanding in two directions 
and contracting in one direction) asymptotic solutions. 
Examples of NSa and NSb solutions are shown in the 
first and the second panels of Fig. 5. The solutions in 
the region marked SI in Fig. 4 lead to singularities where 
A — > (we term this singularity type I) and these sin- 
gular solutions are the same as those appearing in Fig. 2. 
We divided the SI solutions into two classes, Sla and 
Sib. Sla approaches the singularity as pp > 0, qq > 0, 
rr > 0, while Sib as pp < 0, qq < 0, rr < 0. This 
singularity, since it arises because A crosses zero, can be 
overpassed if we introduce a new "time" variable 



dr = dt/A. 



(35) 



Using this r, two solutions Sla and Sib can be joined 
via the singularity, which is shown in the third panel of 
Fig. 5 (Sla,b). The solution Sib, solved backwards in 
time, goes into another singularity which has different 
property from the one between Sla and Sib. At this 
singularity, which we call type II, A goes to -co, and 
at least one of the expansion rates (q in the case of Fig. 
5, Sla,b) diverges. We can say that this solution (Sla 
and Sib joined together) comes regularly from t = — oo, 
turns back at type I singularity, and then goes backwards 
in time into a type II singularity. Or we can also see this 
as two solutions, one coming from t = — oo and the other 
from the type II singularity, "pair-annihilate" at one type 
I singularity. S2 is yet another solution, which comes 
from one type II singularity and disappears into another 
type II singularity. As H avT becomes larger, the boundary 
between Sla and Sib (A = line in Fig. 3) gets pushed 
to approach the line X + Y = 1, and accordingly the 
nonsingular regions NSa and NSb become smaller. This 
is consistent with Fig. 2, which shows the existence of 
the upper limit of H avT for the regular solution through 
the cr < region. 



B. Solutions through the a — 2.5 cross section 

As is expected from the isotropic case discussed in the 
previous section, the solutions through the a > cross 
section are quite different from those through the a < 
section. In Fig. 6 we show the prohibited region (dark 
shaded), A > region (white), and A < region (light 
shaded) in the X-Y plane, with a fixed to 2.5. The el- 
liptic allowed region of small H avr is the one discussed 
in relation to the future asymptotic form of the solution, 
which is expressed as X 2 + Y 2 — XY = 1. As H avl 
takes large values the A = contour takes complicated 
forms, and the region connected to the isotropic solution 
becomes small. 

Figure 7 shows the solutions passing through the X-Y 
plane of the a — 2.5 cross section, and the time evolu- 
tion of each type is shown in Fig. 8. H avl is chosen to 
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be 0.01, 0.05, and 0.1. As H avr increases the nonsingu- 
lar region becomes smaller, and for H avl larger than 0.1 
the nonsingular solution completely disappears from the 
X-Y plane. All the singularities appearing in Fig. 7 are 
type I, and these singular solutions can be extended fur- 
ther by using r defined by Eq. (|35|). Just like the Sla 
and Sib solutions in the previous subsection, Sic and 
Sid, extended beyond the type I singularity, turn back 
futureward and then go into the type II singularity. The 
only difference between these and Sla,b is the direction 
of time, and the former can be regarded as the "pair cre- 
ation" of cosmological solutions, while the latter is the 
"pair annihilation." 

One of the nontrivial results of our analysis, and what 
makes this model very different from ordinary universe 
models, is that the "initial singularity" in the isotropic 
limit is categorized into the type I singularity (see the 
third panel of Fig. 7 and compare Sid and isotropic solu- 
tions in Fig. 8). This means that the singular solutions in 
the model proposed by Antoniadis, Rizos, and Tamvakis 
H or Rizos and Tamvakis Q will, if small anisotropy 
is included, terminate suddenly at finite past with finite 
Hubble parameter or, if extended using r, turn back to- 
wards the future. 

All nonsingular solutions in Figs. 7 and 8 continue 
to the asymptotic solutions expanding in all directions 
in the past asymptotic region. That is, the asymptotic 
form Bl (discussed in Sec. II C) can be reached from 
(7 = 2.5 but B2 cannot. There exist solutions having the 
past asymptotic form B2. For example, solutions through 
the outer white (A > 0) region in the third panel of 
Fig. 3 behave as B2 in the t — > — oo region. These 
solutions, however, go into singularities between a = —10 
and a — 2.5, and do not appear in Fig. 7 or 8. 



C. Nature of the singularity 

In our numerical analysis there appear two types of sin- 
gularities, which we called type I and type II. We discuss 
the nature of these singularities briefly. 

At type I singularities the expansion rates (p, q, r) 
stay finite whereas the time derivatives of them diverge. 
These situations happen when p,q,r ~ \t — t S i ng \ p with 
< p < 1. This is actually the case, and can be verified 
by analyzing solutions near the singularity. Although the 
type I singularity (A = 0) is a physical one, the equations 
of motions can be integrated regularly by using the new 
"time" parameter r, Eq. (|35|), through the type I singu- 
larity. At the singularity p 7 q, r "turn around" (see the 
solutions SI in Figs. 5 and 8), meaning that p, q, and r 
are multiple-valued function of t. Since dt — Adr = at 
the singularity, the solutions are tangential to t = t S i ng . 
Also, in order that the solutions change the chronological 
direction, the leading power of p, q, r in the expansion of 
t near the singularity must be even. Thus we can express 
t using p, q, r as 



t - t„ 



tp,2l{p - Psing) 21 + t p ,2l + l(p ~ p s ing) 2l + 1 + 
tq,2m(q — qsing) 2m + tq,2m+l{q — qsing) 

t r ,2n(r - r sing ) 2n + i r ,2n+i(y - r : 



21+1 

2m+l 
2n+l , 



., (36) 



with l,m,n being positive integers. By solving these with 
respect to p, q, and r, we have 



P = Psing +Pi(t- t smg ) 1/21 + 

Q — Qsing H~ Ql (t ^sing) ^ ~^ 
f ' T 'sing ^1 ~~ ^sing) ^ H~ 



(37) 
(38) 
(39) 



Thus, < l/2Z,l/2m,l/2n < 1, and p = (pi/2l){t 



1/2Z-1 



■, etc., will diverge. The behavior of a 



is similar to that of p, q, and r. In our numerical cal- 
culations (SI of Figs. 5 and 8), l,m,n take the values 
I = m = n = 1, which is the most generic case. 

In the vicinity of other types of singularities, analytic 
expressions of the solutions are obtained by assuming 
power-law behavior of the scale factors and the regularity 
of the modulus field, just as in the isotropic case [[§. 
Because of the anisotropy, there are 3 cases of singular 
solutions other than the type I. 



CI : p ~ pi/t.q^ qo,r ~ r ,a ~ ct ,c7 ~ ai,a ~ er 2 , 
3, 9£ 



Pi = l,<?o + »'o - 77A 

2 oa 



(Jiqor = 0, 



(40) 



C2 : p - pi/t,q~ qi/t,r ~ r ,a ~ a ,a - <7i,a ~ <r 2 , 

<Jir = 0, (41) 



11 3 x ^ 



C3 : p ~ pi/t.q^ qi/t,r - ri/t, a - a ,<r ~ <Jit, <r - <r 1 , 

01 = 0, (42) 



1 d£ 

Pi= qi = ?'i = 1,1 - -A — 

2 00 



where we have chosen the origin of t at the singularity. 
CI is the case where only one of the three expansion rates 
(p, q, r) is singular, C2 two, and C3 all three. The solution 
CI agrees well with our numerical results (S2 of Fig. 5). 
Since Pi,qi,r x = 1 (if not zero), the divergent behavior 
is determined by the sign of t; i.e., if a solution goes 
into a singularity futureward (£ — > —0), then p — > —00, 
etc., and if pastward (t — > +0), p — > +00, etc. Putting 
Eqs. @-@ into Eq. @, we have A -> -00 for 
CI and C2 but A — > +0 for C3. Therefore, according 
to our definition of the singularities (A — > for type I 
and A — > —00 for type II), CI and C2 will be categorized 
into type II and C3 will be categorized into type I. In the 
isotropic case {p = q = r), we can show that A > is 
always satisfied, and C3 can be seen as a type I singularity 
"pushed towards the infinity." 

All of these singularities, both types I and II, are physi- 
cal singularities. This can be shown by putting p, q, r and 
p, q, f into the curvature scalar R — 2(p + p 2 + q + q 2 + 
r + r 2 + pq + qr + rp) . 



G 



D. Summary of numerical results 

We extended the nonsingular universe model proposed 
by Rizos and Tamvakis p| to include anisotropy, and ex- 
amined solutions in the vicinity the of isotropic solution 
in both a < and a > regions. We found both nonsin- 
gular solutions and singular solutions. Nonsingular solu- 
tions inhabit the region where the anisotropy is small, 
and they evolve from the past asymptotic region, su- 
pcrinflate, and then lead cither to Fricdmann-type or to 
Kasner-type solutions in the future. Singularities appear- 
ing in our analysis are classified into two types, namely, 
type I and type II. The type I singularity corresponds 
to the crossing of A = 0, and p, q, and r diverge, while 
p, q, r, and a, a stay finite. At the type II singular- 
ity, on the other hand, A tends to — oo, and p, q, and 
r will diverge. The evolution of the singular solutions 
is characterized by the behavior of A. At the origin of 
the X-Y plane (isotropic solution) A is always positive 
regardless of the values of a or H avl , and as anisotropy 
increases there appear A < re gion s or regions prohib- 
ited by the constraint equation (11), which is shown in 
Figs. 3 and 6. There are three types of singular solutions 
appearing in our analysis (if two branches connected by 
a type I singularity are counted as one solution). The 
first type of singular solution is marked by the crossing 
of the A = line, which we termed the type I singular- 
ity, in the future. This includes Sla and Sib in Figs. 4 
and 5, and if we continue the solution beyond A = by 
changing the variable, this singular solution can be seen 
as a pair annihilation of the A > and A < branches 
of solutions. The A > branch continues from the in- 
finite past, while the A < branch leads to a type II 
singularity at the finite past. The second singular solu- 
tion is very similar to the first one, except it crosses the 
A = line in the past. This solution, examples of which 
are Sic and Sid of Figs. 7 and 8, can be regarded as a 
pair creation of two branches. The singular solution in 
the isotropic model (|^] is a special case of this second 
singular solution. The third singular solution lies always 
in the A < region and never crosses the A = line, 
i.e., includes no type I singularity. This solution is born 
in the type II singularity, and disappears into the type II 
singularity (see S2 of Figs. 4 and 5). 



IV. CONCLUSION 

In this paper we presented anisotropic nonsingular 
cosmological solutions derived from the 1-loop effective 
action of the heterotic string. We found nonsingular 
solutions which evolve from the infinite past asymp- 
totic region, superinfiate, and then continue either to 
Friedmann-type or Kasner-type solutions. The singular 
solutions of moderate anisotropy are classified into three 
types, and involved in these solutions are two types of 
singularities: one corresponds to A = and the other to 



A = — oo. The initial singularity in the isotropic limit is 
a special case of the A = singularity. 

Violation of the energy conditions, which is necessary 
to avoid the singularity, is achieved by the existence of a 
Gauss-Bonnet term coupled to a modulus field. This can 
be confirmed by using the asymptotic forms (|||) and (|3~i| ) 
for nonsingular solutions. We define the effective energy 
density and effective pressure as 



e 
Pi 
Pi 

P3 



-G°o = pq + qr + rp, 
G\ = -(q + r + q 2 + r 2 + qr), 
-(r + p + r 2 + p 2 + rp), 
-{p + q+p 2 +q 2 +pq)- 



G 3 3 



(43) 
(44) 
(45) 
(46) 



Assuming the asymptotic solution A in the region t — > oo 
and that a is always positive (in our numerical analysis 
a keeps its sign except the singular solution S2), the ef- 
fective energy density and pressure behave in the future 
asymptotic region as e, pi ~ iw 2 |i| -2 (i — 1,2,3). Thus, 
the weak energy condition e + pi ~ io 2 \t\~ 2 > and the 
strong energy condition e + J2Pi ~ 2u 2 |i|~ 2 > are 
satisfied. In the past asymptotic region, on the other 
hand, the asymptotic solution becomes B in which the 
Gauss-Bonnet term is dominant. Then e ~ (lu\uj2 + 
w 2 w 3 + L0zuji)\t\~ 4 and p\ ~ — 2(w 2 + cj 3 )|i| -3 — (wf + 
u> 2 +<x>2W3)|i| , etc. If the weak energy condition is sat- 
isfied, e+pi > for i — 1,2, 3, and so 3e+J2Pi ~ —4(u>i + 
L02 +u>s)\t\ must be positive. If the strong energy con- 
dition is satisfied, e + Y^Pi ~ — 4(u>i + w 2 + w 3)|£|~ 3 
must be positive. Neither of them is possible as long as 
H avY ~ (u>i + uj-2 + ^3)|£|~ 2 /3 > 0. Therefore, weak and 
strong energy conditions are violated in the past asymp- 
totic region. 

One of the biggest advantages of our model is that it 
includes a rather long period of (super)inflationary stage 
in a natural form. However, our result, which admits the 
evolution of an almost-isotropic superinflating solution 
into a Kasner-type anisotropic solution (see the panel 
"NSb" of Fig. 5 for example), suggests that the superin- 
flation in our model does not isotropize the space-time. 
Our recent study 0] of the cosmological perturbation for 
the homogeneous and isotropic background shows the ex- 
istence of exponential growth in graviton-mode pertur- 
bation during the superinflationary stage. Together with 
this we can conclude that the cosmic no-hair hypothesis 
H does not hold in this superinflationary model. We are 
interested in whether this result is common to all kinetic- 
driven superinflation. 
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-0.5 1 1 ' 1 ' 1 ' 1 1 

-10 -8 -6 -4 -2 2 4 

FIG. 1. The H-a phase diagram of isotropic solutions 
(A = 1, H > 0, a > 0). Time flows from left to right since 
a > 0. The nonsingular solutions are plotted with a solid 
line, and singular solutions with a dotted line. The bold line 
is a critical solution marking the border of singular and non- 
singular solutions. All solution flows in the H > 0, a < 
quarter-plane continue smoothly to the H > 0, a > quar- 
ter-plane. In the H > 0, a > quarter-plane, however, only 
the flows below the critical solution continue to the a < 
region. 
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FIG. 2. (a) The average expansion rate in the anisotropic 
case. The equations are solved from a — —10, where the 
initial anisotropy is fixed as X = 0.1, Y — 0.2. For large 
initial H avr , there appear singularities with which the solution 
flows terminate suddenly, keeping H avr and a finite, (b) The 
behavior of -Havr, -Havr, o, a, and a in a singular solution 
appearing in (a). Initial conditions are the same as in (a) 
except the initial -Havr is set to 0.04. We can see -H aV r and a 
diverge, but H avl , a, and a stay finite. 
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ffavr = 0.005, a = -10 
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ff avr = 0.01, o-= -10 
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FIG. 3. Constraint and A on the a — —10 section of 
the X-Y plane. H avr is 0.001, 0.005, 0.01 from above. 
A > 0, A < 0, and excluded regions are indicated by white, 
light-shaded, and dark-shaded areas, respectively. In the dark 
gray region the Hamiltonian constraint ( |ll| ) is not satisfied, 
and A cannot be defined. Cosmological solutions inhabit the 
white and light gray regions, and those in each region are 
separated by sin gula rities since p, q , r become infinite when 
A = [see Eqs. @, @, and @]. 
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H avr = 0.001, a = -10 



H avl = 0.005, a = -10 



# avr = 0.01, a = -10 






0.2 0.4 0.6 0.8 1.0 1.2 X 0.2 0.4 0.6 0.8 1.0 1.2 X 0.2 0.4 0.6 0.8 1.0 1.2 X 

FIG. 4. Solutions through the a = —10 cross section, ffavr is 0.001, 0.005, 0.01 from the left. NS means nonsingular. NSa 
leads to a Friedmann-type solution (expanding in all directions) and NSb leads to a Kasner-type solution (expanding in two 
directions and shrinking in one direction) in the future asymptotic region. SI means it leads to a singularity where A — > 0. Sla 
is the solution whose behavior near such a singularity is pp > 0, qq > 0, rr > 0, while Sib behaves as pp < 0, qq < 0, rr < 0, 
near the singularity. S2 leads to a singularity where A — ► — oo. 
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FIG. 5. Behavior of solutions appearing in Fig. 4. t = is the time when a — —10. These are solutions through the points 
(X,Y) = (0.1,0.2), (0.2,0.4), (0.4,0.6), and (0.8,0.9), respectively, in the H avr = 0.01, a = -10 plane. 
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ff avr = 0.1, (7 = 2.5 




X 
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H avl = 0.2, cr = 2.5 
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H avI = 0.3, cr = 2.5 
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FIG. 6. Constraint and A on the a = 2.5 section of the 
X-Y plane. H av r is 0.1, 0.2, 0.3 from above. A > 0, A < 0, 
and excluded regions are indicated by white, light-shaded, 
and dark-shaded areas, respectively. 
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FIG. 7. Solutions through the a = 2.5 cross section, -Havr is 0.01, 0.05, 0.1 from the left. NS means nonsingular, which leads 
to an expanding universe in the past asymptotic region. SI means it leads to a singularity where A —> 0. Sic is the solution 
whose behavior near such a singularity is p > 0, q > 0, r < 0, while Sid behaves as p > 0, q > 0, r > 0, near the singularity. 
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FIG. 8. Behavior of solutions appearing in Fig. 7. t — is the time when a — 2.5. The initial values are chosen on the 
a = 2.5 plane as (H avT ,X,Y) = (0.01,0.6,0.8), (0.1,0.6,0.8), (0.1,0.1,0.2), and (0.1,0.0,0.0) for NS, Sic, Sid, and isotropic 
examples of the solution, respectively. 
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